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Abstract 

We investigate the role of the transverse spin dependence in Drell Yan lepton 
pair production to NLO in QCD at parton level. In our analysis we deal with 
the large px distributions. We give very compact expressions for the virtual 
0{a 2 s ) corrections to the cross section and show that the singularities factorize. 
The study is performed in the MS scheme in Dimensional Regularization, and 
with the t'Hooft-Veltman prescription for 75. A discussion of the structure 
of the real emissions is included, and detailed methods for the study of these 
contributions are formulated. 
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1 Introduction 



There has been considerable recent interest in the study of polarized DIS collisions 
at HERA, both from the theoretical and the experimental side. Most of the work, 
so far, has been directed toward the analysis of the (longitudinal) spin content of the 
proton, in connection with the so called "spin crisis" and the puzzling results of the 
EMC collaboration (see 0,0)- The next step of this ongoing analysis is to study, 
for various ranges of x and Q 2 , the leading twist (polarized) parton distributions in 
p-p collisions in a different experimental setting. Such a program will be possible in 
the not distant future at polarized hadron colliders such as RHIC, using polarized p-p 
incoming states. 

A direct gluon coupling, absent in DIS, will allow us to study the gluon content 
of the nucleons in an independent manner. The coupling is non-anomalous and, 
therefore, all the theoretical ambiguities associated with the scale-independence of 
the axial anomaly are not present. 

There is also much more to be gained from these studies. 

It is well known that in polarized DIS, the leading twist transverse spin distribution 
hi(x) decouples, being chirally odd, and one possible way to get some information 
about its parton content is to resort to polarized p-p collisions. 




Figure 1: The Drell-Yan process 

There are at least two reasons to justify a detailed NLO study of the Drell Yan 
(DY) process, Fig. [I], both in the case of the longitudinal and of the transverse asym- 
metries. The first one is that the final state is particularly "clean" , due to the presence 
of a lepton pair. This means that the photon fragmentation contribution does not af- 
fect the momentum distribution of the final state. In particular, the scale dependence 
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associated to a fragmenting photon is absent. 

The second reason is that transverse asymmetries are sensitive to gluon distribu- 
tions only through radiative corrections. In previous works, the complete (real plus 
virtual) NLO corrections to the non singlet DY cross section with longitudinally po- 
larized initial states have been presented (for the nonzero pt = qr distributions). In 
this work we make the first step toward extending these calculations to the case of 
transversely polarized incoming states. 

Although much of the experimental data are usually collected at small p?, the 
study of the distributions at large pt remains, nevertheless, a challenging and complex 
problem. 

As a results of these studies, we are able to perform nontrivial perturbative checks 
of factorization in QCD and we developed more complete strategies for the investiga- 
tion of spin dependent cross sections. Most of the techniques that we consider in this 
and in former related || work are, in fact, applicable to the easier case of polarized 
e + e~ collisions. 

Our analysis is organized as follows. After a discussion of the main features of 
the various cross sections in the the study of the Drell Yan process, we analyze the 
properties of the cross section da/d 4 q dfl. 

We separate the contributions to this cross section into "diagonal" and "non- 
diagonal" terms, and then we give very compact results for the virtual corrections 
to the "diagonal" terms. The term "diagonal" refers to those contributions obtained 
by contracting the hadronic tensor (at parton level) with the g^ u part of the leptonic 
tensor, while the "non-diagonal" part contains the momenta of the lepton pair (here 
denoted as v%v%). We show that these "diagonal" virtual corrections factorize. The 
result is presented in a remarkably simple form. The reason for performing this 
separation comes from the structure of the factorized cross section in terms of parton 
distributions as given originally by Ralston and Soper in the parton model @. We 
will come back to this point briefly in one of the following sections. 

Then we move to a detailed discussion of the real emission contributions and show 
that in our cross section the integration over 3 of the five final-state partons is sufficient 
to cancel the singularities of the virtual corrections. We will not give here simplified 
results for the real emissions, which will be presented elsewhere, but we set up the 
stage for a complete analysis of the real emissions to NLO. Some relevant technical 
aspects of the calculation have been included in various appendices. A very efficient 
method to perform the renormalization of all the diagrams -which is of simple symbolic 
implementation- is illustrated in some detail, since it is of general applicability. In a 
concluding section we summarize our analysis and elaborate on future developments 
and unsolved issues of the theory. 
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Figure 2: Leading order contribution to the Drell-Yan process 



2 Zero distributions. General Features 



Let Hi and H 2 be two polarized hadron that collide thereby producing a heavy photon 
and some hadronic remnants (X) 



H l + H 2 ^ 1 * + X. 
At parton level the process goes through the elementary scattering 

a + b -»• 7* -y l + 



(1) 



(2) 



and the leading channel is the annihilation channel. The process is purely electro- 
magnetic at the lowest order. In fact, Fig. |2| is the leading order contribution to 
the DY process. At this order, the lepton pair is produced with a total zero qx, i.e. 
back-to-back. 

If we include radiative corrections -with an additional parton in the final state 
beside the photon- then the lepton pair gets a nonzero qx- Two cross sections are 
studied in general: cr(Q 2 ) and da/dydxp, with y being the rapidity of the photon and 
xp is its "Feynman's x" . Both cross sections are sufficiently inclusive and well defined 
in perturbation theory so that an analytic calculation is possible. 

The NLO corrections to the cr(Q 2 ) cross section are quite straightforward to obtain. 
At this order, for instance, the real emissions involve only 2 partons in the final state, 
while the virtual require the 1-loop renormalized q q 7* vertex plus quark self-energy 
insertions (see Fig. 




Li 



'COCCCCOT 



— www 

Li 



Figure 3: a s order contributions 

In this cross section one of the final state partons is integrated over (the gluon 
or the quark) and the transverse momentum dependence of the lepton pair is not 
completely resolved. 
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A better accuracy is obtained by going to NNLO, i.e. at order a 2 . 

For instance, in the non-singlet sector, all the diagrams Vi,...,Vn (Fig. 7) are 
needed, plus there are contributions from the so called "G", "F" and "H" sets (see 
Figs. 4-6). Basically, the G-set involves real emissions of the form qq — > gg'y*; the set 
F involves both annihilation (qq — > qq')*) and scattering (qq — > qq^y*) diagrams, while 
the set H involves qq' — > qq'j* scattering diagrams. 

Additionally, one needs to evaluate all the corrections coming from the interference 
between the 2-loop (polarized) on-shell quark form factor and the lowest order process 
of Fig [| The latter are the only two-loop contributions needed for a complete NNLO 
calculation of the invariant mass distribution. 

In a completely anticommuting scheme for 75 this can be obtained directly from 
the literature, and differs from the result of refs. |T2j just by an overall helicity factor 
(1 — h). In the t'Hooft-Veltman scheme it is expected that additional helicity- violating 
terms will appear, since the on-shell form factor is not infrared (IR) safe. However, a 
redefinition of the polarized splitting function of the form AP qq = P qq + ACF(-(z — 1) + 
0(e 2 ) should be sufficient to eliminate the helicity- violating contributions. However, 
as we have mentioned before, the invariant mass distribution has to be analyzed both 
in its diagonal and non-diagonal contributions in order to obtain information on the 
perturbative corrections to the transverse spin cross section. In fact, as we will discuss 
below, it is the structure of the parton model result for a TT that forces us to consider 
a more general approach to the DY process, and renders the calculation far more 
complex than in the unpolarized or in the longitudinally polarized case. 



3 Nonzero pt Distributions 

Since there is a lot of information in the transverse momentum distribution of the 
lepton pair and since, experimentally, one needs to set cuts below a certain p T , the 
study of the transverse momentum distributions turns out to be of considerable inter- 
est. Various cross sections are studied for this purpose. One of the most studied cross 
section is defined as da/d 2 QT, or, equivalently as da/dt du. In the analysis of the vir- 
tual corrections we will concentrate on the latter. We remind that the diagonal part 
of the more general cross section da/d 4 q dQ, with Q defined to be the solid angle in 
the rest frame of the photon, is easily extracted from da/dt du. In their non-diagonal 
contributions the two cross sections are, of course, different. 

The distribution cross section da/dt du starts at order a s and the two diagrams L\ 
and Z/2 denote its two Born level contributions in the non-singlet sector, and proceeds 
to NLO (or 0(a 2 s )) through the radiative corrections of Figs. |J-[7|. 

The virtual and the real corrections to this cross sections are the same corrections 
that appear in the NNLO analysis of a(Q 2 ), the invariant mass distribution. 

The extension of da/dtdu to the small qx region can be performed by a resum- 



mation [20| in impact parameter. At large qx, instead, the cross section is perfectly 



5 



well defined. 

If we use Dimensional Regularization (DR) (with n = 4 — 2e) to calculate both 
the real and the virtual corrections to this cross section, then the overlapping soft and 
collinear regions of integration will produce double poles 1/e 2 and the collinear re- 
gions single 1/e poles. The study of this cross sections in the non-singlet sector shows 
that "factorization" takes place both for the unpolarized 0] and for the longitudinally 
polarized case ||. By "factorization" here we mean that the double poles will cancel 
between the real and the virtual contributions, while the left over collinear singu- 
larities "factorize", i.e. can be absorbed into universal parton distributions, which 
are process independent. The usual factorization picture involves also a collinear ex- 
pansion of the momenta of the hard scatterings around a light-cone direction, fixed 
by the large momenta of the two colliding hadrons, and arrested at a certain order. 
Parton distributions of various twist, in this picture, are defined as matrix elements 
of nonlocal operators at light-like separations. Natural ways to incorporate into this 
picture contributions of higher twist and proofs of factorization have also been 
formulated. 

An important part of the analysis is the proof of the gauge invariance of the 
factorization formula up to a certain twist. The operatorial expansion of shares 
the same features of ref . |T3 , although the analysis is performed directly on the light- 



cone, rather than in momentum space [19[ 



In both formulations, the violations to the scaling behaviour of the twist 2 (and 
of any twist) structure functions are purely logarithmic, at any order in perturbation 
theory and all the hard scatterings -for distributions of any twist- are calculated in 
the collinear approximation. 

These approaches, which are manifestly gauge invariant, do not isolate any trans- 
verse (explicit) qx in the dominant cut-diagram describing the factorization formula. 
Isolation of a fixed qx in the parton distributions brings the operatorial expressions 
that describe these distributions away from the light cone. 

In the following, we will assume that the leading cut-diagram describing the Drell 
Yan process at nonzero pr contains a hard scattering of the form L\ or L2. Therefore, 
in the leading twist approximation, the qq partons undergo hard scattering with purely 
collinear light-like momenta. The qr of the photon, therefore, is not "intrinsic", as 
discussed in [0, |i~If , but it is due to the real (lowest order) emission of a gluon. The 
scaling (In qj) violations to this basic picture are generated by the additional radiative 
corrections of Figs. 4-7. 



4 cr(Q 2 ) versus da/dtdu 

In order to obtain cr(Q 2 ) from da/dt du, then we need to integrate the latter cross 
section in t and u. By doing so, we will generate additional singularities (up to 
1/e 4 ) after integration, since the transverse momentum distribution contains matrix 
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elements of the form 1/t, 1/u, 1/ut etc, which will be responsible for higher order 
poles in e. These higher order singularities, according to factorization, will cancel 
in the final sum of real and virtual, after adding the contributions coming from the 
two-loop on-shell quark form factor. 

Therefore, the 0(e°) radiative corrections to the cross section da/dt du gives, 
after t — u integration, only the bulk of the total NNLO corrections to the invariant 
mass distribution. These corrections, in the longitudinally polarized case, have been 
calculated ||. 



5 da/d A q dQ and transverse spin at nonzero Pt 

In their pioneering paper on the polarized Drell Yan process, Ralston and Soper 
introduced the notion of transverse parton distributions. It is by now well established 
that hi(x) is a leading twist (or twist 2) distributions which can be measured in p-p 
collisions using transversely polarized beams. 

As we are going to see, in order to have a non vanishing projection on the transverse 
spin distribution of the hadronic tensor, the parton model prediction requires that the 
lepton pair should not be completely integrated out. The result is common to all the 
approaches 0, [| HH followed so far in the analysis of the parton distributions in 
polarized collisions. 

Here we briefly overview these results in order to make our treatment self con- 
tained. 

The factorized (parton model) picture of the collision emerges in an infinite mo- 
mentum frame, with a large boost parameter (formally P — > oo). 

In this limit the two colliding hadrons have large "plus" or "minus" components 
and fields are effectively quantized at equal light-cone time x + . The degrees of freedom 
of a fermion are halved, with the negative ip_ = l/27~7 + (bad or dependent) light- 
cone projections of the Dirac fermions functions of the positive ip + = 1/27+7^ (good 
or independent) ones ||| [IIJ, and of the transverse components of the gauge field. The 



plus components can be labeled by helicity eigenvalues and chirality eigenvalues 

of the same sign. The negative projections are basically interpreted as quark-gluon 
composites. The requirement that the helicity of the composite be ±1/2 forces the 
composite to be of the form if)]/ 2 ~ %p + 1 ^ 2 A +1 and ip^ 1 ^ 2 ~ ip^^A' 1 . The chirality 
of the bad quantum components have flipped sign compared to the helicities. While 
the distributions fi(x) and g±(x) have a natural interpretation (i.e. are diagonal) in 
the chirality/helicity basis, the distribution h\(x) does not. Since a fermion of fixed 
transverse spin is in an eigenstate of the Pauli-Lubanski operator, it is convenient 
to expand the components of the fermion in a basis of eigenstates of such operator 
(transversity basis), which are expressible as linear combinations of the helicities. 
hi(x) then measures the number of quarks polarized in the transverse direction minus 
the number of those polarized in the opposite direction in a transverse polarized 
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hadron (at a fixed light-cone momentum fraction x). 
One starts by defining the hadronic tensor 



= J ^^(PaSa; PbS b \ [J m (0), J v {x)) \P A S A - P B S B ), (3) 

with P A and Pb being the momenta of the colliding hadrons and Sa and Sb their 
covariant spins. In a collinear basis 
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2P B 



S A = -^^n + - -^n 

Sa = '^UL n - - l^l n + , (I) 
M 2P B 

with spin vectors orthogonal to the hadron momenta. 

Pa ands P B are the two large light-cone components of the incoming hadrons, of 
equal mass M. 

Assuming factorization, the quantum numbers of the two hadrons can be decoupled 
by Fierz transformations, and the interaction described at leading order by a single 
hard scattering, with a single photon in the unitarity diagram. 

The distribution functions that emerge -at leading order- from this factorized 
picture are correlation functions of non-local operators in configuration space. They 
are the quark-quark and the quark- ant iquark correlators. 

Their expression simplifies in the axial gauge, in which the gauge link is removed 
by the gauge condition. For instance, the quark-quark correlator takes the form 

(® a/A ) ap (P,S,k) = J ^e ik -*(P,S\4\Q)^\z)\PS). (5) 

We have included the quark flavour index a and an index A for the hadron, as usual. 
Fields are not time ordered since they can be described by the good light cone com- 
ponents and At, as discussed in 0. 

Using hermicity, parity and time reversal invariance, one can give a general ex- 
pansion for these correlators in terms of basic amplitudes which are functions of P,S, 
k and are Dirac algebra valued. 

However, the basic matrix elements that appear in leading power factorization are 
not the correlators themselves, but their integrals over the k~ and /cr components of 
momentum. 
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Specifically, one introduces the Sudakov expansion 



k 2 - k 2 

k = xp + n + H —n~ + kr (6) 

2x 



and may decide to integrate over both k and kx in (||). 
Using the identity 

d(k-p)d 2 k T ik . z r dX 



e lk ' z = / ^-S A (z - \n)e lXx , (7) 



(2vr) 4 J 2ir 

and projecting over the Dirac basis of the 16 independent 7 matrices, one obtains the 
general expansion of the correlator 

J(f> = l~jTr[$\- 275 / Tr + Tai - j Tr [7"$] 

+ 757 ^Tr [ 7At75 $] + i l5 a^ J Tr [z^ 75 $] , (8) 

where the sign of integration / is a shorthand notation for the k~ , kx integration in k 

J Tr [T$] = - j d(k ■ p)d 2 k T Tr [T$(k)] 
r d\ _ 

= J ^(p,s\mm^)\ps). (9) 

Therefore, after integration, the correlators are sampled on light-cone rays. 
The transverse spin distribution h±(x) appears in the projector of the quark-quark 
correlator over 20^75 

47T 

+MA(/n 1 ' - p u n»)h L (x) + M 2 (S T n u - S^)h 3 (x). (10) 



Using Jaffe's definition of twist [10] t = 2+number of M, one finds that h\(x) is 
twist— 2, hi{x) is twist— 3 ands h^{x) is twist— 4. 

The leading twist expansion of the quark-quark correlator then is of the form 

^e iX *(PS\m^(\n)\PS) = ± {i>h{x) + A75 fai(x) + 75 fir t>h x {x)) . (11) 

In their parton model analysis the authors of || obtain expressions for the hadron 
tensor W^ u in terms of parton distributions which are peaked at qx = 0, (eqs (67), 
(68) and (69) of §). 
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A more general approach has been discussed in [IT]. There, the authors, following 
Ralston and Soper, assume a factorization with an intrinsic qx for the total momentum 
of the partons entering hard scattering. The factorized expression for the hadronic 
tensor is assumed of the form J7| [ll| 



W" = \Y. 6 b- a e 2 a J d 4 k a d*k b S\k a + k b -q) 

ab 



xTr[<5> a/A (P A S A , K)-f® h/B {P Bi S B , k b )Y}- (12) 

Notice that this primordial is not induced by radiative corrections. In this last 
approach, by setting q = (or q T = 0) one is yet unable to recover the results of Jaffe 
and Ji, since the integrals over the parton momenta k A and k B in the hard scattering 
are still coupled. They can be uncoupled only if one starts with a general qr in the 
factorization formula Eq. (|i~2l), and integrates over q?. Therefore, the approaches of 
[|TT| , [|] compare at the level of the integrated cross section da/dq + dq~dfl, an expression 
derived originally in 0. The authors of [[11] also claim to be in disagreement with 



the expression quoted in M for the Alt or longitudinal-transverse asymmetry. 

The parton level results for the g^-integrated cross sections, however, in both 
approaches, agree, and are also in agreement with the result given by Ralston and 
Soper, that we simply quote 



x (l-cos 2 9)(l- \ A X B gi(x A )g 1 (x B )) 

+ sin 2 8 cos(20 — <fi A — <fi B ) \S A \ \S^\hi(x A )hi(x B ). (13) 

The angles 9 and in Eq. (p~3|) refer to the leptons and are measured respect to 
the Collins-Soper (see below) frame, while <p A and <p B are the azimuthal angles of the 
two transverse spins 4- vectors in the same frame. 

It clearly shows that an integration over the azimuthal angle of the lepton pair 
makes the projection over h\{x) vanish. 



da 

dq + dq~dfl 



6 Parton Level Analysis 

Since the structure of the final state in the analysis of the transverse spin cannot be 
neglected, it is convenient to set up some diagrammatic notations and divide all the 
calculations into two subsets. Diagrammatically, this is easily understood in terms 
of cut graphs. For instance, the contraction of the hadronic tensor with the diagonal 
part (i.e., with g^ u ) of the leptonic tensor 

L^u ~ v lfi v 2 u + v 2lJ ,v lv - g^Vx ■ v 2 (14) 
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is described in terms of cut diagrams of cross sections for the process qq —>■ 7* . 

The non diagonal part, instead, is schematically described by cut diagrams of cross 
sections for the process qq — ► l + l~ . 

In the next sections we will analyze the NLO structure of the non-diagonal non- 
singlet virtual corrections for the px distributions. 

Notice that if we integrate completely over the final leptons and use the gauge 
invariance of the partonic matrix elements, then only the diagonal projection survives. 

In our calculation we use Dimensional Regularization (DR) all the way. The chiral 
projectors are regulated according to the t'Hooft-Veltman prescription (HV) ||15| . 

There are clearly many advantages in being able to use this regularization, given 
the complexity of the calculations. However this is not generally possible in most of 
cases, except for specifically defined cross sections. As we have just mentioned, the 
cross section has to be sufficiently inclusive in order to be able to cancel the infrared 
divergences order by order in perturbation theory. At the same time, in the case of 
transverse spin, we need to leave in part the lepton pair in the final state unintegrated. 
As we are going to show in the next sections, one possible choice is given by the cross 
section da/d 4 q dfl, with Q evaluated in the rest frame of the photon, which we are 
going to investigate. 

7 About the treatment of 75 

Here are some remarks concerning the extension used for 75 in n dimensions. In the 
t'Hooft-Veltman prescription for 75 for n = 4 — 2e dimensions, the Dirac algebra is 
split into a 4 and into ann-4 subspace. 

The metric space is also split in a similar way 

9^ = dpu + 9^, (15) 
and the Dirac algebra written direct sum 

7m = 7 M + %• (16) 
The anticommutation relations for the 7 matrices are 

The antisymmetric tensor e^ upa is 4 dimensional and projects only over the 4 dimen- 
sional subspace 

^pvpaQaf] 0, tfii/pa'Ya 0. (1^) 

75 is defined in n-dimensions as 



[%>%] 



2g 



0. 



(17) 
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75 = ^e^YYYY- (19) 

With this definition, 75 anticommutes with the 4-dimensional 7's and commutes 
with the remaining the n-4-dimensional ones. 



8 Distributions at nonzero 



We calculate the spin dependence of the cross section using the relation [14" 



u{p,s)u(p, s) 



-p 1 (1 - A75 + 75^±) 



(20) 



-I. 
1. 



where s± represents the transverse components of the spin vector s M , and s^s^ 
Clearly, A and |sj_| are not independent, and in fact satisfy the relation A 2 + \ s±\ 

In our NLO calculation a s — > a s ([i 2 ) is the running coupling constant renormalized 
at the scale /x 2 in the MS scheme, satisfying the renormalization group equation 



with 



,da s 

A) 
ft 



-a. 



11JV, 



G 



34iV 2 



Air J 



) 2 + 0(a 



2N P 



10N C N F 



(21) 



- 2C F N F 



(22) 



and with iVV the number of active quark flavors. 

The lowest order (non-singlet) contributions to the real emissions arise from the 
two, q + q — > 7* + g, amplitudes L\ and L2. As mentioned before, they are also the 
lowest order contributions to the pr distributions in Drell Yan. 

We refer to these two diagrams as the direct and the exchange (or crossed) ampli- 
tudes, respectively. 

The squares of the direct amplitude L\ and exchange amplitude L 2 in Fig. § in 
n = 4 — 2e dimensions are given by 



M dd (sx,s 2 ) 



M cc (s 1 ,s 2 ) 



2 2 2 Cp 2u 

e f 9 g s 



Nr. t 



h{l 



2s 



-e" ( si ± ■ s 2 ± + —(k 2 ■ s 1± )(k 2 ■ s 2± ] 
tu 



000 C p 2t 

e f 9 9s 



N r . U 



[l-e) 2 -h(l + e) 
2s 



-(? I Si ± ■ S 2± + —(k 2 ■ S 1± )(k 2 ■ S 2± ) 

tu 



(23) 
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e/ the charge of the quark, and we have set A1A2 = h. The quantity Cf/N c is the 
color factor, and a s = g 2 s /Ait is the QCD strong coupling constant. 
The interference term is more complicated and is given by, 



2M dc ( Sl ,s 2 ) = 

r ■ [(1 - e)(Q 2 s - etu) - h(l + e)(Q 2 s + etu) - 2hetu 
Nc tu L 

+(s 1± ■ s 2L )((Q 2 s + tu) - (1 + e)(Q 2 s + etu)) 
{fa ■ s 1± )(fa ■ S 2 J_) 



2 2 2 4 

e f g 9 



+ 



Q 2 



(2sQ 2 (l-e-e 2 )) 



The sum of the two Born amplitudes squared is 

\M B ( Sl , s 2 )\ 2 = M dd ( Sl , s 2 ) + 2M dc ( Sl , s 2 ) + M cc { Sl , s 2 ) 

= z)?&%Tu I (1 ~ m ~ e) ( 2g2fi + (1 " e)( ' 2 + u2) ~ 2etu ) 

+4he {Q 2 s + (u + t) 2 ) 

+(s 1± ■ s 2 ±) (2tu - 2e(s + u)(s + *)) - e 2 {t + u) 2 ) 

The distribution cross section is defined by 
da B 



dt 



= ^ 



K 2 is defined by 



C 



Ti(Q 2 , u, t, h) + T 2 (Q 2 , u, t, s ± i, s ±2 ) 

1 fA7Tfi 2 Y fsQ 2 Y 



tu 



(24) 



(25) 



(26) 



(27) 



'N c T(l-e) \Q 2 ) 

where we have rescaled, a s — * a s (/i 2 ) e , so that it remains dimensionless in n = 4 — 2e 
dimensions. The first component of the cross section is the same as longitudinally 
polarized cross section 



Ti(Q 2 ,M,t, h) = (l-h)T B (Q 2 ,u,t) + 8he 
where T B is the unpolarized born level matrix element. 

T B (Q 2 , u, t) = 2(1 — e) 



Q 2 s +(u + tf 
tu 



^ e) |! + i) + W-'i 



t U 



ut 



2e 



(28) 



(29) 
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The transverse spin contribution to the cross section is 

2tu - 2e(s + u)(s + t)- e 2 (t + u) 2 



T 2 (Q ,u,t } s 1±} s 2± ) = 2s 1± -s 2± '- 



+4(fci • s 1± )(h ■ s 2X ) 



tu 

s{2{l-e)tu-e 2 {t + u) 2 ) 



t 2 u 2 

tu / tu 



W { 2(1 - e) - e 2 ^±^) - Ae(s 1± • s 2± )^-, (30) 



where 

2(fci ■ sij_)(fci • s 2 ±)s 

P = : \-si±-s 2 ±. (31) 

tu 

If Si± = s 2 ± = 0, Eq. (|26|) reduced to the longitudinally polarized cross section cross 
section of Ref . || . 
In the limit h = 

• s 2± = —si± ■ s 2 ± = cos (32) 

and 

2{h ■ s 1± )(h ■ s 2 ±)s 2(kt ■ s 1± )(ki ■ s 2± )s 

" h Si_L • S 2± = Si ± • S 2± 

tu tu 

2k 2 s 

= -^cos6»icos(6»i + AO) -cosA6» 
tu 

= cos(2#i + A#). (33) 

Here, k\ = (tu)/s and 6\ is an angle between k± and s*i_l. 

Away from 4 dimensions Ti clearly contains helicity violating contributions, due 
to the prescription chosen for 75. 



9 The Structure of the Non Singlet 

The radiative corrections to the non singlet come from the interference between the 
two diagrams L\ and L 2 and the diagrams denoted Vi, ....V\i 

^virtual _ 2R( , + L J £ y*\ ? ^ 

The other real qjqf diagrams are 

^t-lE^ + lE^I 2 , (35) 
1=1 1=1 
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where the amplitudes Fi are shown in Fig. 4 and the amplitudes Gi are shown in 
Fig. 5 with / = /'. The absolute square of the F; amplitudes can be written as: 



i=i 



E^l 2 = 

Fx + F 2 \ 2 + \F 3 + F 4 \ 2 + 2Re (]T F t £ F*) + | £ F t \\ (36) 



\i=l i=5 / i=5 

having used 

Re((F 1 + F 2 )(F 3 + F 4 )*) = 0. (37) 
Therefore, the diagonal real g/g/ contributions are 

da r g e f g l f = d&i + da 2 + da 3 + da F , (38) 

where 

da, ~ iFi + F^ + I^TG,! 2 , (39) 

i=l 

d<r 2 ~ |F 3 + F 4 | 2 , (40) 



da" 



3 



2i?e(£F i ]TF*V (41) 



cia F ~ lE^I 2 - ( 42 ) 



i=5 



The real diagonal g/ qf diagrams are 

^SS~IE^-E^I 2 . (43) 

i=l i=4 

where the amplitudes Hi are shown in Fig. 6 and the relative minus sign between the 
direct and exchange diagrams is due to Fermi statistics since / = /'. The absolute 
square of the Hi amplitudes can be written as follows: 

4 8 

IE^-E^I 2 = 

i=l i=4 

| E Htf + I E ^l 2 - 2i?e ^ E AT) • (44) 

i=l i=4 \i=l i=4 / 

Thus the diagonal real g/g/ contributions are, 

da r q e f a q l f = da n + da 4 , (45) 
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where 

1 



da H ~ lE^r + lE^ ( 46 ) 

i=l i=4 



da 4 ~ -2ite 

Vi=l i=4 / 



(47) 



The diagonal (/ = /') real part of the non-singlet cross section is given by subtracting 
1~3D from fl38T) as follows: 



d&NS = + + da 3 - da 4 , (48) 
where we have used the fact that dap = dan which arises due to 



/ \i:m 2 = U ix>i 2 +i/ ix>r- (49) 

JPS 3 ^ 2 JpS 3 ^ 2 J PS, ^ 



When integrating over the two identical particles £3 and over the phase space PS3 
in ( j64|) an extra statistical factor of 1/2 must be inserted so as not to double count. 
Thus, the diagonal (/ = /') part of the complete non-singlet cross section is given by 

d&NS = dai + d&2 + da 3 — d&4, (50) 

where 

da x = da( act + da virtua \ (51) 
and where da( act is the "factorized" cross section given by 

da[ act = dai + da{ ounter . (52) 

It is necessary to remove the collinear singularities from the initial state by subtracting 
the cross section d^ unter to da\. 

The off-diagonal real q/qp diagrams are 

^t,~iE^i 2 , (53) 
1=1 

where the amplitudes F{ are shown in Fig. 4 and where / 7^ /'. Similarly, the off- 
diagonal real qpqp diagrams are 

^,HE^I 2 > (54) 

i=l 

where the amplitudes Hi are shown in Fig. 6. In this case, 

IE^I 2 = IE^I 2 , (55) 
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and these two contributions cancel when forming the non-singlet cross section. 

The total contribution to the transversely polarized cross section can then be 
written in the form 

da = da NS + 2 (da H + da 4 ) ff + + da r q e f f fi ) , (56) 

where the suffix // refers the diagonal (in flavour) contribution and the remaining 
term in bracket are the off diagonal scattering diagrams in the F and H sets. 

10 Off-shell renormalization and tensor reduc- 
tions 

Although in the non-singlet sector we do not encounter any anomalous diagram, the 
regularization scheme still suffers from unphysical features, such as helicity violation, 
which are absent in other schemes. We should also mention that in the case of the 2-to- 
2 contributions, all the dependence on the hat-momenta of the the hard scattering can 
be eliminated. This is equivalent to ask that the momenta which are not integrated 
-and also the spin 4-vectors- have only 4 dimensional components. 

The calculation of the radiative corrections to the vertices is often performed by 
using various tricks in order to simplify the loop momenta, such as the Feynman 
parameterization combined with symmetric integration. An alternative way to pro- 
ceed is to use the Passarino-Veltman (PV) |18| recursion procedure in order to relate 
the tensor integrals which appear in the calculation to the corresponding scalar ones. 
The PV method introduces scaless integrals which have to be handled with particular 
care. In an appendix we illustrate in a simple way how to obtain the renormalized 
expression of all the coefficients of the expansion of tensor integrals to scalar form. 
The approach || is easy to implement in symbolic manipulations. We start from the 
expression of any tensor integral and apply the Passarino Veltman (PV) reduction 
procedure to determine the coefficients of the expansion. Only the renormalization 
of the scalar self energy diagram B (Q 2 ) (and of its Q 2 — > limit) is needed. We 
use dimensional regularization with a single parameter e to regulate both UV and IR 
singularities. Notice that scaless integrals such as B (0), so generated, are intrinsi- 
cally ambiguous. While it is possible to set such terms to zero by definition (on-shell), 
off-shell renormalization of B(Q 2 ) leaves us with a 1/e pole contribution in B (0) after 
that the Q 2 — > limit is performed. We isolate and remove the UV singularity (by 
setting n = 4 — 2e), we then switch e — ► — e to regulate the remaining IR divergences. 
Finally, we send Q 2 — * 0. This leaves us with a 1/e pole for Bq(0), of IR origin. 
Details of the method can be found in the appendix. 
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11 NLO Diagonal Virtual Corrections 



In this section we present results for the diagonal corrections to the NLO non-singlet 
cross section and show their factorization. Although the evaluation of these diagrams 
is a formidable task, the result can be condensed in a remarkably simple form. We 

get 



sda 



virtual 



a. 



dtdu 



e}K 2 ^-5(s + t + u - Q 2 ) (T x + T 2 ) 



1 - 



a r(l -e) /4vr/i 2 
27rr(l-2e) \Q 2 ~ 



+ ^-[(l + h)R 1 + R 2 ]Y 



where 



3tu 



-2C F 

+ [l 

- \2Li 2 



/8(2Q 2 s + t 2 + u 2 



Q 4 s(t + u) 



tu(t + u) 
t 2 + u 2 



\ tu 
(2C F - N c )Li 2 (- 

/ ' ^ • l.r 



t«(s + t)(s + w) (s + t)(s + M) 

* + u\ 2Q 2 s + t 2 + u 2 + 2s 2 



t-Q 2 



' t-Q 2 



N, 



tu 

2s + t 



c- 



+ 2d 



s 2 + {s + uf 



u 



tu 



+{2C F -N c )[2\n{-^ ) U ^ A2 

- In 2 



g 4 - (t + m) 2 + ln2 r_8_\ s 2 



(t + u) 2 \Q 2 J tu 

f\t\\\ s 2 + (s + u) 2 ' 



tu 



AQ 2 s - 2st + tu 
{s + u) 2 



+ N ( 



t 



c 



s + u 



at ^ (ItW, (\ U W 2Q 2 s + t 2 + u 2 \ u 



(57) 



(5? 



describe the NLO dependence on the longitudinal components of spin. The trans- 
verse spin contributions are condensed into the expression 



R2 = (Sl± • S 2 ±) 



» A r ^ +2(20 ,_ Wc) ( ln m,_ 



2s , / \u 
+ ln 



2s 



Q 2 s + t / 



{AC F + N C )^-2C F (§ + 



2s 



t + u/ \ 6 tu 
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+ {2(2C F -N c )Li 2 ' ' + 'M-s-'+'r 



s J u 2 



■2 



t-Q 2 \ \ 2C F (s 2 + u 2 ) - N c (s 2 - u 2 




Q 2 J V V t u s + u 



(59) 



Notice that in the rest frame of the photon, (3 = si± ■ s 2 ±, and the expression above 
simplifies substantially. From Eq. ( |57D it appears obvious that the virtual corrections 
factorize. In fact the double poles 1/e 2 correctly multiply the 2-to-2 Born cross section 
calculated in n dimensions, evaluated with the most general spin dependence in the 
initial state. If we send any of the transverse spin vectors to zero, then the result 
coincides with the longitudinally polarized corrections obtained in ref. ||. Finally, 
by sending to zero any of the helicities and any of the 2 spin vectors, one reobtains 
the result given in ref. in the unpolarized case. Eq. ( p7[) clearly contains (finite) 
helicity violating contributions which are not canceled by the real emissions in the 
MS scheme. Finite subtractions are needed in order to restore helicity conservation 
at parton level [HJ . A discussion of the longitudinal contributions to the real emissions, 
which are a subset of the total emissions considered here -and briefly analyzed below- 
can be found in [HI. 



12 The qg sector and the Asymmetries 

When the incoming quark and the incoming gluon are both transversely polarized, 
the scattering matrix element is zero. 

On the contrary of the longitudinal case (LL scattering), in which the process 
takes contributions from the non-singlet and from the singlet sector, TT scattering in 
Drell Yan involves only the quark sector. 

This statement can be directly verified to lowest order by an inspection of the 
Dirac traces, and remains true to all orders in perturbation theory. 

To be definite, we introduce 2 (transverse) polarization vectors for the gluon, here 
denoted as e±± and e 2 j_, and choose a collinear base made out of the 2 momenta pi 
and &2, the momenta of the initial state quark and of the final state quark or gluon. 
P2 is the gluon momentum. All these partons are massless: pi = p 2 , = k\ = 0. 
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The two transverse polarizations four-vectors are then defined as 



1 



e 2 ±n = i = ==€ VM favy$kX. (60) 

\jl{k\ ■ h)(pi ■ fa)(pi ■ fa) 

One easily gets 

h±. = ( A /&2 i>2- A^2 • P2+ Mpi ■ P2- fopi ■ fa)^ 5 (61) 

which can be easily implemented in symbolic calculations. With these definitions then 
e± • fa — e± • p 2 — e± • fa — 0, which guarantees transversality of the (physical) gluon. 

For initial state asymmetries, in order to have a nonvanishing trace, each spinor 
projector has to contribute with a zero or an even number of gamma matrices. The 
contribution ~ 75 jb is clearly odd and therefore renders the total number of gamma 
matrices odd (vanishing) in each trace. 

A longitudinal quark, however, can scatter from a transverse gluon, since only the 
~ 75A part of the fermion projector becomes relevant. These contributions are part 
of the LL scattering cross section, since we can relate transverse gluon polarizations 
to the usual gluon helicities by a linear combination e ± = l/v / 2(e 1 j_ ± ie 2 ±)- The 
vanishing of Att in the qg sector indicates that the ratio between All and Att is 
1 to lowest order. This simple result shows that Drell Yan plays a key role in the 
analysis of the transverse spin distributions, since the gluon contributions is absent 
to all orders. 



13 The treatment of the real emissions 

In this section we focus our attention on the analysis of the real emissions. The 
steps that we describe here are necessary in order to proceed with an analytical or a 
numerical study of the real diagrams. 

There are 2 basic modifications to be considered in this case: 1) the presence of 
a 2-to-4 final phase space integration; 2) the presence of transverse spin in the HV 
regularization. Some of the features appear already in the longitudinal case, but we 
prefer to present a comprehensive discussion of all these aspects for future reference. 

Let v\ and v 2 denote the momenta of the 2 leptons in the final state. We also set 
k-y — q for notational convenience and introduce the invariants 

s = ( Pl +p 2 ) 2 , 

512 = (h + k 2 ) 2 , s 23 = (k 2 + k 3 ) 2 , 
h = (Pi ~ h) 2 , Mi = (p 2 - h) 2 , 

513 = (fa + k 3 ) 2 , 
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h = (pi - h) 2 , h = (pi - h) 2 , 
u 2 = (P2 - h) 2 , u 3 = (p 2 - h) 2 , 

t' i = ( Pl ~v i ) 2 , u'i = (P2-Vi) 2 . (62) 
Only five of them are independent since 

s = Q 2 - ( Ul + u 2 + u 3 ) = Q 2 - (ti + t 2 + t 3 ) = sia + s 23 + sis - 2Q 2 , 

S12 = S + U 3 + t 3 , 

S23 = S + Mi + ti - 2Q 2 , 

s 13 = s + u 2 + t 2 , 

h = Q 2 + t\ + t' 2) 

Ul = Q 2 + u \ +u ' 2 . (63) 

The integration of the 2-to-4 matrix elements introduces -beside the usual matrix 
elements of the LL cross section- new matrix elements containing t\, t' 2 , u\ and u' 2 . 

In the separation of the 2-to-4 into 2-to-2 subintegrals, we choose to work in two 
separate frames 1) the center-of-mass frame of the k 2 — k 3 pair; 2) the center-of-mass 
frame of the photon with axis fixed by the Collins-Soper [^TJ choice. 



To be specific, we start from the 4-particle phase space integral is given by the 
formula 

PS 4 = J ^ 1 ^25 + (^)^+(^)^ n fc3^fc25(A;3 2 )5(A; 2 2 )5 n (p 1 + p 2 - k 3 - k 2 - v x - w 2 ).(64) 
We lump together the momenta k 3 and k 2 in k 23 and V\ and v 2 in q = k± as follows 



PS 4 = J d n qd n k 23 6 n ( Pl +p 2 - q - k 23 ) J d n Vl 5 + {vl)^q ~ vif) 

x / d n k 2 6 + (k 2 2 )6 + ((k 23 -k 2 ) 2 ). (65) 



From this expression it is obvious that in the rest frame of the (2, 3) pair (k 23 rest 
frame) the collinear singularities linked to the two gluons can be easily isolated in 
Dimensional Regularization. 

If hat-momenta are not present, then we easily get 



PS,= d n q d n v l 8 + {v 2 1 )8+{{q-v 1 ) 2 ) d n k 23 5 n ( Pl +p 2 -q-k 23 ) 



x 



W 2^r[n/2-3/2] X h2 (66) 

I 12 = r d9i sin n " 3 Q x d6 2 sin n " 4 6 2 (67) 
J o J o 



where 
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is evaluated in the rest frame of k 2 z and from which the singularities can be exposed. 
The final result can be cast in the form 

PS, = J d n q ^ d n Vl 5 + (vf)S + ((q - ViY) 4 2 - 2 2n _ 2r>/2 _ 3/2] x I 1>2 , (68) 

where we have used the relation s + t + u = S23 + Q 2 ■ The remaining left-over 
integral has to be evaluated in the photon rest frame, with the Collins-Soper choice 
of axis. Finally (|68| ) can be expressed in the form 



d n q dn n n -\Q 2 ) n ' 2 - 2 s% 2 - 2 ^j— — x h 2 . (69) 

y ) 23 2"- 2 r[n/2-3/2] ' v ; 

Matrix elements containing t\ , such as t\/(t 2 u 2 ) and similar, can all be exactly 
evaluated by this procedure, slightly generalized to the tensor case. Matrix elements 
containing k 2 ' ^1 (and any denominator) are also easily handled in these two frames. 
After integration, the result can be conveniently expressed in the Collins-Soper (CS) 
frame. For this purpose [^TJ we parameterize the momentum of the photon in a 
collinear basis defined by the two incoming partons p\ and p 2 , and expand q = Xipi + 
X2P2 + <?t- We introduce two metric projectors and g q defined by 



i 1 v v A 1 

Jiv to/ PlP2 PlP2 



9t =9 



Pi -P2 Pi' P2 



< = sT - (70) 

These two metrics project in the directions orthogonal to the collinear basis and 
to q respectively. 

An orthonormal set of 4-vectors (q,z,x,y) is then defined as follows 
M Z ix X fi y» 

y Q 4^ V^X 2 V-Y* 



where 



\Pl-P2j \Pl-P2j 
\Pl-P2j \Pl-P2j 

Y" = -^—e^p lv p 2p q a , 

Pi -P2 

Pi, = - ^V- (72) 
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The set is orthogonal and spacelike. In the photon rest frame, therefore, it defines 
a cartesian set of axis. Note that in the collinear frame of the two incoming parton, 
Z is collinear to p\ and p 2 , while projects over the transverse (w.r.t. pi and P2) 
plane. 

The vector V\ can be expanded as 



< = y + -| (sin6>cos0x M + sin 9 sin </> y* 1 + cos0 z^) (73) 

and the cross section expressed in a form closer to eq. ( |T3"D . Finally, one expresses the 
scalar products p\ • v 1, q ■ V\ or s±i ■ V\ in the CS frame. 



14 Hat-momenta Integration 

As we have mentioned before, the HV prescription introduces matrix elements con- 
taining hat-momenta, due to separation of the n-dimensional space into 4 and n — 4 
dimensions. Beside v±, the only other hat-momenta expected from the matrix ele- 
ments are those of k 2 and k 3 , which are treated in the way discussed below. Other 
hat-momenta can be eliminated by a convenient choice of axis. This causes, at least 
in part, a modification of the phase space integral which appear in the unpolarized 
case. 

In the case of unpolarized scattering, the singularities are generated by poles in the 

matrix elements which have the form l/t 3 , 1/W3, 1/(^3^3) and similar ones, in multiple 

combinations of them. Multiple poles can be reduced to sums of combinations of 

double poles by using simple identities among all the invariants and by the repeated 

use of partial fractioning. This is by now a well established procedure. In our case we 

encounter new terms of the form l/t\ and I/M3 and new matrix elements containing 

typical factors of the form fc 3 , k 2 , and k 2 ■ k 3 at the numerator. In the longitudinal case, 

the hat-momenta that appears in the matrix elements are those related to pi,p 2 ,k\ 

and k 2 . As we discuss in the Appendix, we can set to zero all the matrix elements 

containing hat-momenta of pi,p 2 and k\ by a convenient choice of axis. Since the 

separation of the n-dimensional space into a 4 and an n-4 dimensional subspace is 

arbitrary, we can always assume that all the momenta which are not integrated over 

are embedded in the 4-dimensional part. Therefore, in the case of a general spin 

dependent cross section, we can remove all the hat-momenta related to pi,p 2 , and 

s 2 ±. The left-over hat-momenta, which are not set to vanish, are those related to 

Vi, k 2 and k 3 , while v 2 is reexpressed in terms of q and Vi. Since q is not integrated 

over, q is also zero. There are 2 ways to integrate the hat-momenta contributions. 

The most direct one is illustrated below. For instance, in the integration of k 2 , we 

2 ^ 

use the on shell condition k\ = to relate k 2 to the 4-dimensional projection k 2 . 
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We start from the matrix element in the (2, 3) rest frame 

PS 2 = j d n k 2 5(k 2 2 )5((k 23 - k 2 ) 2 )t 2 . (74) 

We have set k 2 = (k 2 , k 2 ), with 

k 2 — k 2 (l, cos #3 sin 6 2 sin 6 1 , cos 9 2 sin 0!,cos #i) (75) 
being the 4- dimensional part of k 2 . We easily get 

k 2 = Sj f sin 2 #3 sin 2 6 1 sin 2 2 . (76) 
Therefore the usual angular integration measure 

n-2 

dn {n - 2) = [] sin""'" 2 9 t dO t (77) 
1=1 

with 

Q (n-2) = 2 n T sin™-'- 2 ^ ^ (78) 

r=i Jo 

is effectively modified to 

3 n-2 

dfi (n " 2) = JJ sin"-' 9 1 d9 l x J] sin"-'- 2 0,d0,. (79) 
1=1 1=4 

The intermediate steps of the evaluation are similar to those in the previous section. 
In n = 4 — 2e dimensions we get 



PS, = 



ir 2t e fut- Q 2 s 23 



2% 4 r[l - e] 



[ 7T d6 l sin 6l~ 2t I* d6 2 sin6 2 2 ~ 2t , (80) 
2 Jo jo 



where 0\ and $2 are the only relevant angles which appear in the matrix elements 
and therefore are not integrated. We have displayed also the #3 integral since it is 
different from the unpolarized case. 

The contributions from hat-momenta are either finite or of order 1/e and therefore 
play a key role in the cancelation of the mass singularities in the cross section. 

15 Collinear Subtractions for da^s 

The real emission diagrams (sets G,H,F) have collinear singularities which are gener- 
ated by the emission of a massless parton off the initial state quark or gluon. Factor- 
ization in the quark sector involves diagrams from all the 3 sets mentioned above. As 
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discussed in j| for the unpolarized case and in g in the case of longitudinal polar- 
izations, the diagrams which require explicit factorization, in the non singlet sector, 
are those of the set G. The remaining collinear singularities cancel in the other real 
sets after adding all the contributions. For this purpose we start defining the relation 
between the bare and the renormalized structure functions by the equation 



G A ^(x,M 2 
with R of the form 



1 dz_ 

z 



a. 



In 



ribare 



R 



l An - Jl-e 
M ; T 1 -2e 



'47T/i 2 



+ C^-(z) 



51) 



(82) 



for the LL subtractions and 



*A P (z) T[1 - e] 



'Ann 2 



+ <V ff (z) 



53) 



for the TT (or transverse spin) subtractions, where the finite pieces C^j(z) are 
arbitrary. M is the factorization scale. 

In the MS scheme they are set to be zero. The longitudinal, AP^, and the 
transverse, A T P gq , splitting functions are given by 



AP, 



qq 



P 



qq 



C7 



(1 + z 2 ' 
(1 - z). 



-6(1 
2 v 



AP, 



mi 



AP 



99 



N, 



c 



'1 + z* 



+ 



33 - 2N F r/1 
+ —^5(l-z) 



AP gq = C F [2-z], 



A T P qq 



2x 



K (l-x)+ 
with the distribution 1/(1 — z 



1 . f(z 
dz- 



defined by 

1 /(«)-/(!) 



l-z 



15) 



Beside the usual subtractions of the singularities in the longitudinal spin sector, 
now we have also similar subtractions in the transverse spin sector. 

As we have mentioned above, the collinear subtractions, in da^s involve only the 
set G. In these diagrams a quark (or an antiquark) can emit a collinear gluon from 
the initial state. R^j is therefore proportional to AP qq (= AP qq ) or to AtPjq. 
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The non singlet cross section is given to order a 2 by the renormalized cross section 
d&i with the the collinear initial and final state singularities subtracted 

daf act da x 
s— — - — = s 



dtdu dtdu 

a r l da^ LL 

/ dz l R q ^ q {z 1 ,M 2 )s f- \ Pl ^ ZlPl 5(z 1 (s + t-Q 2 ) + u i 

ATT Jo at v ' 

a r 1 d<j( 1)LL 

^E/ dz 2 R q ^ q {z 2 ,M 2 )s^^- \ P2 ^ Z2P2 5(z 2 (s + u-Q 2 )+t) 



- 7T dz x R rq ^ q {z x , M 2 )s — f- — \ pl -+ zlP i 6(zi(s + t-Q 
zn Jo at v 



u 



~ T^E/ dz 2 R q ^ q (z 2 ,M 2 )s^f— \ P2 ^ 2P2 5(z 2 (s + u-Q 2 )+t%86) 
In k Jo at 

In Eq. ® da LL and da TT are Born level cross sections at a rescaled value of one 
of the two incoming momenta ( Z\P\ or Z\p 2 ) due to the collinear gluon emission, and 
evaluated in n dimensions. 

LL and TT refer to their longitudinal/transverse spin content and are proportional 
to Ti and T 2 respectively. 



16 Conclusions 

We have presented a general discussion of the spin dependence of the Drell Yan cross 
section for the nonzero pt distributions, and shown that the analysis of the hard 
scattering can be performed in the t'Hooft-Veltman scheme consistently. Factorized 
expressions for the virtual corrections have been presented. Our analysis has been 
focused on the non singlet sector, which is the main source of transverse spin depen- 
dence in this process. The result for the virtual corrections presented here, however, 
are those of the entire process. 

From our results, we find that -at non zero pr- a transverse spin dependence is 
already generated from the diagonal part of the leptonic tensor. The calculation is 
an explicit proof of factorization of the hard scattering , in a nontrivial case. We 
have, along the way, developed a complete methodology for the analysis of the real 
emissions in the same scheme. 

We should also mention that factorization theorems for spin-dependent cross sec- 
tions are yet to be proved. In fact, all the results presented in the literature are 
obtained assuming a light-cone dominance of the process. Particularly important - 
and unsettled- appears the issue of higher twist effects in the factorization formula 
and in the expression of the Alt asymmetries. 

Work toward a complete analysis of the hard scattering for the nonzero pt distri- 
butions is now in progress and we hope to return on some of these issues in a future 
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work. 
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17 Appendix. The 2-particle phase space 

Let's start considering the virtual contributions to the cross section. The relevant 
phase space integral is given by 

It is convenient to introduce light-cone variables (q + , q~ , q±), with 

q = q + n + + q~n~ + q± (88) 

and 

n^-^O^il) (89) 

and work in the frame 

pt=P2=\[\ (90) 

to get 

PS* = J^^\ql\^-M q l\8(2 q + q --\ q l\-Q 2 )8(s + t + u-Q^ (91) 
= J d + qd-q(2q + q-Q 2 r^ 2 6(s +t + u- Q 2 )j^Z~ 2 (92) 

with 

f7r n-3 ry nl2-\ 







using 



and after covariantization 



d(q + ,q ) = J_ 
d(t,u) 2s'' 



(94) 



(2q+q- - (?) -+ (Q 2 -tKQ 2 -)-*Q 2 = * (95) 
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PS 2 = Jdtdu-{-) r(n/2 _ 1} ^ , (96) 

* = T^T / 1^4^ f i£Y f *( a + t + u - Q 2 )UQ\ u, t). (97) 



AN C J (2s) 2 \Q 2 J \ut 
The integration to the total cross section can be obtained as follows. We start 
from 

1 f s -Q^ n - 3 

W 2 - 1 2 2 "- 4 r[n/2 - 1] 
We also define 

J= I d9sm6 n - 3 d9. (99) 



/ dPS 2 = .„ - - * r , ( s n ' 2 - 2 r (sin 6) n - 3 d9. (98) 

J 7r n/2-l 2 2n-4r 7l/2 - 1 I S / Jo^ J V ' 



We define as usual 

s = (Pi + P2) 2 i = (pi - ^i) 2 u = (p 2 - h) 2 k\ = Q 2 (100) 
with s + t + u = Q 2 and Q 2 is the invariant mass of the photon 

We define the Bjorken variable x = Q 2 /s (0 < x < 1) and we assume that x is 
not parametrically small (x ~ 0) or close to 1, x ~ 1. Both the small-x region and 
the large-x (x — > 1) need resummation, since large logs in the variables x and 1 — x 
are generated. The scattering angle is cos 9 = (Q 2 — 2u — s)/(s — Q 2 ) and u varies in 
between — (s — Q 2 ) < u < 0, (i. e. — s(l — x) < u < 0). 

It is convenient to introduce the angular variable y = 1/2(1 + cos#) and rewrite 
)8|) as 

l PS > = 2t^tt f^T" W1 " »>i" /M - (101) 



We get 



J = 2"- 3 / 1 dy(y(l- 2/ )r/ 2 - 2 . (102) 

In terms pf the variable y we have it = — Q 2 /x(l— x)y and t = —Q 2 /x(l—x)(l—y). 
We set n = 4 + 2e and define 

Using these notations -PS^, ] = s]J[, ]. 
We get 

x 2 T[n/2 - 2] 2 



PS 2 [l/(ut)] = K[n, s] 



(Q 2 ) 2 (l-x) 2 I>-4] 



P9 ri/fl - P9 \l/v] - ~ X K\n d r[n/2-2]r[n/2-l] 
PS 2 [u/t] = PS 2 [t/u] = K[n, s] T[n/ l~ 2 ™ 2 \ (104) 
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which can be expanded up to the desired order in e. For instance, the real emissions 
from the qq — > 7* sector give 



a 



'2ue 2 e 2 f 



1 - e 



2N, 



c 



a s C F \ /r(l-e) 



(1-h) 



'4(5(1 -x) 4 1+x 2 



e(l-x) 



Q 2 ) V 4tt J \Y(1 -2e) 
+ 8(l + x 2 ) 



V (1-x) 



41n(x) 



1+x 2 
(1-x). 



+ I6/1 



-6(1 -x) (5(1 - x) 1 - x + x T 



2e 



, / -4(5(1 - x) 2(1 - x) 2(3 - 2x + 3x 2 ) 
+ (si± • s 2 _l) 1 h 



-45(1 -x) +2(l-x)ln 



(1-x). 



(1-x) 2 



(105) 



We have factorized the expression in this manner to facilitate comparison with the lit- 
erature || , in which the normalization convention is to remove the first two quantities 
in parentheses in Eq. (|105| ). 



18 Appendix. Integration of the real emission di- 
agrams 

In order to integrate over the matrix elements, we need to evaluate the various scalar 
products which appear in such matrix elements, in the cm. frame of the pair (1,2). 
For this purpose we define the functions 

A(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2yz — 2xz 
\ l > 2 (x,y,z) 



P[x,y,z] 



2^ 



E[x,y,z] = X+ 2 V v - Z . (106) 



It is easy to show that 

\px\ = P[s 23 ,pt,ui 



2 1 S23 - Ui 



2^ 

P2I = P[s23,Pl,tl 



.2 + 1 _ g 23 ~ ^1 

2^ 



\fo\ = \fa\ = P[S23,PI,P 2 ] 



2 ' 
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We define 



We get 



P°i 



pI 



\k\ = ^P[s^] J K ^\ . 

V ^23 2^/s^ 
Pi = ^[s 2 3,P?,Ml 

P° 2 = E[s 23 ,plt 1 ] = S -^^- (107) 



S23 


- Ml 




/S23 




-h 


2 V 


/ S23 


t - 


Q 2 



s 

tt, = -7T7h?- (108) 



S 2 3 - Ul sv 



s 2 3-h s(l-v)-Q 2 



Q 2 - s + S23 s(l-v + vw)-Q 2 



'S 2 3 V S 23 

< <■ -4^- (109) 



In the derivation of ( |109| ) we have used the relation s+ti+ui = Q 2 +s 2 3. There are 
four different parameterizations of the integration momenta which we will be using. 
In the first one, which is suitable for unpolarized scattering one defines (in the cm. 
frame of the pair (1, 2)) in general 



! '1, cos 6 2 sin 61, cos#i v 



2 



^2 = 23 (1, — cos 9 2 sin 0x, — cos0i), 

Pi = p\{\, 0, 0, sin^i, cosV>i), 
p 2 = p!j(l,0, ...O,sinV'i,cos'0i), 

ki = (hi, 0, ...0, \k±\ sin ^2, \ki\ cos ip 2 ), (HO) 



where the dots denote the remaining n — 2 polar components. We get 

sin^j 



s(l-w) \ 1/2 



1 \s(l-vw)-Q 2 J ' 
. smip(Q 2 - s(l - v -vw)) 

3m * 2= — m^r^> — ■ 
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It is convenient to use the parameterizations 
• set 1 



set 2 



set 3 



Pi = p?(l,0,..., 0,0,1), 

p 2 = p°(l, 0, sin 0,008^"), 

k\ = (fe°, 0, —\ki\ sinip, 0, |fei| cos?/'), (H2) 



Pi 



Pi(l, 0, sinip" , 0, cos-?/'")) 
p 2 =p°(l,0,.., 0,0,1), 

fei = 0, I fei| sin^', 0, |fei| cosi/j'), (H3) 



Pi = p?(l,0, ...,sinV>,0,cosV>), 

P2 = °> ■■•> - sin ^ 0> COS ^'), 

fei = (fe?,0, ...,0,0,1^1), (114) 

where 0, ... refers to n — 5 components identically zero. It is straightforward to obtain 
the relations 



cos ip" 
cos ip = 
cos ip' - 



_ (-S23 - h)(s 2 3 ~ Ul) ~ 

(S23-*1)(S23-Ul) 

(-g 2 + s - g 23 )(g 23 - ^i) - 2g 23 (g 2 - h) 

X 1 / 2 (s,Q 2 ,s 23 )(s 23 -u 1 ) 
(s 2 3 - ti)(-Q 2 + -s 23 ) ~ 2s 23 (Q 2 - Ul ) 
X 1 / 2 (s, Q 2 , s 23 )(s 23 - ti) 



(115) 



We select a specific set depending upon the form of the hat-momenta in the matrix 
elements. We choose the center of mass of the gluon pair (2, 3). Following the notation 



of reference 22, 231 we define 



4 k ' l) = r dd x r d 6 2 sin (a + b cos 6 X 
Jo Jo 

x (A + B cos 6»i + C sin 9\ cos 2 )~ l , 



(116) 



where a,b, A, B, C, are functions of the external kinematic variables . 

This expression trivially contains collinear singularities if k > 1, I > 1 and a 2 = b 2 
and/or A 2 = B 2 + C 2 . The singularities are traced back to the emission of massless 
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gluons from the initial state and in the final states. In this special case, it is convenient 
to rescale the integral and let the angular variables ip, ip', ip" defined above appear. 
The cases a 2 = b 2 and A 2 ^ B 2 + C 2 , a 2 ^ b 2 and A 2 = B 2 + C 2 and a 2 ^ b 2 and 
A 2 7^ B 2 + C 2 are discussed in ref. p2| . It is easy to figure out from the structure 
of the matrix elements which cases require a four dimensional integration and which, 
instead, have to be evaluated in n dimensions. This procedure is standard lore. Notice 
that only two independent angular variables 6\ and 6 2 appear at the time. 

Matrix elements containing more than 2 integration invariants (for instance 
\jt 2 u 2 u 3 or l/t 2 t 3 u 2 ) have to be partial fractioned by using the Mandelstam rela- 
tions. Just to quote an example, we can reduce these ratios in a form suitable for 
integration in the 9i, 9 2 variables by partial fractioning 



(117) 



1 1 / u 2 + u 3 \ 

t 2 u 2 u 3 t 2 u 2 u 3 \Q 2 - s - ui) 

where we have used the identity 

Q z — s — u\ 

If we let /[.] denote the corresponding angular integral, we get 

I[l/hu 2 u 3 ) = — — (I[l/t 2 u 2 ) + I[l/t 2 u 3 ]) . (119) 

Q z — s — U\ 

The integrals on the rhs of this equation are now in the standard 8\, 9 2 form. 
Defining 

t] . /' r , Wi r°' 39l T: 4 '' 2 — Trnm- (120) 

Jo Jo (1 — cos#i)*(l =f cos x cos b 1 ! sin x cos ^2 smViy 



we get 



V n [- -] = C[n, i,j]F[l,j, n/2-1, cos 2 1], 

V n [+, -} = C[n, 2, j]F[l,j, n/2 - 1, 1/2 + sin |], 
V n [-+]=V n [--], 

V n [+,+] = C[n,iJ]F[lJ,n/2-l,sm 2 Z], (121) 



where 



l[n — 2 — 1— j\ i A \n/2 — 1J 

Use of the sets 1,2,3 and 4 allows us to get rid of a large part of the hat-momenta 
involved in the real emissions. As an example of integration over the hat-momenta, 
we consider the matrix element 
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= / d n v x 5 + {vl)5( + {q - Vl ) 2 ) [ d n k 2 5 + (k 2 2 )5 + ((k 23 - k 2 ) 2 )^^, (123) 
J J t 2 u 2 

and the sub-integral in the (2,3) frame can be expanded covariantly 

J t 2 U 2 

The evaluation of the coefficients A B and C can be carried out by standard 
methods. However, since we project the n-4 dimensional components by the metric 
g al3 and since pi = p 2 = q = 0, then we easily find that the contribution of this matrix 
element is zero. Therefore, the only hat-momenta left over are those involving k\ and 
v\ (i.e. k\ and v\) which can be treated as discussed in section 13. 



k 2 ■ v 2 
hu 2 



19 Appendix. A fast way to renormalize 

In this appendix we illustrate a simple method developed by two of us || to obtain 
renormalized expressions for any Feynman diagram at one loop. The partons are 
set on-shell from the beginning and UV and IR divergences are identified. After 
performing the PV reduction, a large number of massless tadpoles beside self-energies, 
scalar triangles diagrams and scalar box diagrams, are generated. We identify a set of 
prescriptions for handling massless tadpoles and implement them symbolically. They 
can all be derived by renormalizing off-shell, i.e. by subtracting the UV divergence 
while keeping the external lines off mass shell. We omit the derivation of these rules, 
which have been applied before || and just describe here some simple applications. 
The renormalization of the other scalar diagrams is performed as usual, by subtracting 
the UV poles whenever they appear. 

19.1 Self-energy diagram 

Let's start from the scalar self-energy contribution which is given by 

d n l 1 



hip) 



2 



(2n) n P(l + p) 

-i{^)- nt2 {- V 2 Y 12 i v 1 J , (125) 



u r(2 + cj 

where n = 4 + 2u. The PV function Bo(p 2 ) is defined by 



B (P 2 ) = ^^h(p). (126) 



in 2 
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In the PV reduction, isolated B (0) terms are renormalized off-shell, to give 

Bo-BnO) = -> ( 127 ) 

UJ 

while terms containing the product uj x B(0) are set to vanish 

ujB (0) = 0. (128) 

For off-shell self-energy diagrams we proceed exactly in the same way We renormalize 
isolated self energy contributions which then take the form 

~ 2-\n(^Pj + 0(uj) (129) 

but we leave unrenormalized all the expressions containing the self energy times an uj 
factor 

ujB (p 2 ) = -1+0(uj). (130) 

The latter result can be obtained by using the unrenormalized expression of the self- 
energy 

' /i z J V uj, 

instead of the renormalized one. We will see that these prescription allow us to 
reproduce all the one loop vertices needed in the calculation quite straightforwardly. 
It is also easy to check that the coefficients of the reduction of the (tensor) box 
diagrams are finite, as they should, although the PV reduction introduces spurious 
singularities. 

19.2 Scalar triangular vertices 

The general scalar vertex integral is 

and PV function is defined as 

Co(pI,pIpI) = ^f^2(pi,p 2 ), (133) 

where pz = Pi — P2- 



Bo™^) = ( -£) (2--) (131) 
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In our paper, we need two different types of PV scalar vertex functions. 
• Case (1) p\ = p\ = and p\ = q 2 + 0. 

Pi is the momentum of the incoming quark, p 3 = pi + p 2 is the momentum of the 
gluon and p 2 is the momentum of the virtual quark. 

= (13 4) 

If we set p\ = q 2 the PV function can be reduced to 



q 2 { fi 2 J [to 2 12 / 



1 L 2 6Z& - 7T 2 

I — — H - — 

q 2 V cu 2 u 12 



2 I 2+^ + ^7— ), (135) 



where L g2 = In (-^r) ■ 

• Case (2) p\ = Q,p 2 = Q 2 and pi = q 2 . 



Pi the momentum of the incoming quark, p 2 is the momentum of the photon and 
Pz is the momentum of the virtual quark. In this case 



The PV function is 



= <?r?l^r^ + ^ir^)- < 137 > 

19.3 gluon(photon)-quark-quark vertex (type 1) 

As an application of the renormalization procedure we illustrate here the derivation 
of the vertices which are needed in the calculation. Notice that if we use the method 
discussed above we don't have to use any symmetric integration in order to isolate 
the coefficients of the scalar expansion. 



J d l l 2 (l + m) 2 (l + P,) 2 - (138) 



ITT 2 J l 2 (l+Pl) 2 (l+p 2 ) 

By tensor reduction, we can set 

1 dn h 2 (l+p%+ P2 ) 2 = 9 a "Coo+P a lP fC u +p^ 2 C 22 + (p?pj +p?rf)C 12 . 

(139) 
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We also introduce the expansion 

(27T/i) 4 - n 



> / dni mTT wn ^2 = Pi Ci + P2C2. (140) 

27T 2 7 / 2 (/ +Pl) 2 (l +P2) 

1) We take p\ = 0,j% = Q 2 and (pi -p 2 ) 2 = 0. Then Eq. fllggp becomes, 

V" = (2 -n) {((2 -n)C o-# 22 -^ 2 AC 12 ) 7 M + 2^C 22 + 2^C 12 } 

+(2 - n) (2^ - p 2 27 ^) C 2 - 2 (2hp1-hhY) C x . (141) 

Notice that after applying the PV reduction we will encounter terms of the form Bq(0) 
and ooB (0) which must be handled with the rules given above. In fact we get for the 
tensor coefficients 



^ = 2^-2^-0,(0,0,^), 
B (0) B (Q 2 ) 

2 " q 2 Q 2 ' 

_i 5 (g 2 ) 

L'OO — t + 



4 4 ' 



C„ = ^-3^ + 3^ + C„(0,0,Q 
B„(0) B„(Q 2 ) 



C*22 
C*12 



2Q 2 2Q 2 ' 
_L_ gp(0) £ (Q 2 ) 
2Q 2 2Q 2 2Q 2 ' 



(142) 



The vertex correction becomes, 



W = (2 - n) ((2 - n)C 00 - Q 2 C 22 - Q 2 C 12 ) Y + [n - 2)Q 2 ^C 2 + 2Q 2 ^C X 



2 2L Q2 - 4 



7T 



2 



= 7 T" - 8 + - + 3L ,-L^j. (143) 

The vertex correction in Fig. |3] is this type. 

2) Similarly, let's consider the case of p\ = p\ = and t — (pi — p 2 ) 2 = p 2 . 
In this case, 

fl (0) , B (t) 

l ~ — ~ + ~r ' 

_ fio(O) , 5 (t) 

C2 - — + ~r> 
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r 1 , Bo(t) 

L'OO — T + 



Cxi 

C*22 
C\2 



4 4 ' 

^o(O) B (t) 

2t 2t 

B (0) B (t) 
2t 2t 

1 

2t 



(144) 



The vertex correction is 



V» = (2-n){((2-n)C 00 + tC 12 )r + 2fop%C 22 + 2^C 12 } 

+ (2 - n)2fap£C 2 - 2 (2forf + t^) C x . (145) 

We get 

V = r (£_4 + L«) + ^ pf (l-5 + 2L t ) + ^K(i-l + i. 

= 7 ,(L 4+i( )_^ K (i_ 2+it ) + ^(r_ 1+it ), (146) 

where L t = ln(— t/fi 2 ). Vertices of diagrams V\ and Ves in Fig. |7| are in this category. 
3) v\ = vl = Q 2 and t = (p x - p 2 ) 2 = p§. Pi • p 2 = (Q 2 - *)/2 and 
^>2 + 2pi • P2 In this case the coefficients are 

go(o) w ^o(Q 2 ) , , ,x gog) n2 Co(o,t,g 2 ) 

Ul ~~ 77^2 IT — ^ 77n2 772 + W + *v77T2 772 ~~ °» 



2 



(g 2 -t) - (g 2 -t) 2 ^ '{cp-tf - (g 2 -t) 
g (g 2 ) £o(t) 

(Q2 _ t) ( Q2 _ f) » 



_ 1 , ^ B (Q 2 ) B (t) 

^00 — T + V 77775 77 — r " 



4 " 4(g 2 -t) 4(g 2 -t)' 

Cn = ^ m >_ t) B^ ^B^)) 



(g 2 -t) 2 v y 2(g 2 -t) 2 (g 2 -t) 
-(sg 2 + 4 g 2 t - t \ * o( *> + g 2Co(0 ' *' g2) 



3 



2(g 2 -t) 3 " (q 2 -*) 1 



_ g (g 2 ) B (t) 

^22 — ~7TF7^ — + 



2(g 2 -t) 2(g 2 -t)' 
r 1 n 2 b (q 2 ) 2 B (t) 

12 ~ 2(g 2 -t) ^ 2(g 2 -t) 2 + ^ 2(g 2 -t) 2 ' 1 j 

The vertex correction becomes 
V = (2 - n) { ((2 - n)C 00 - Q 2 C 22 - (Q 2 - t)C l2 ) ^ + 2{> 2 f 2 C 22 + 2^C 12 } 
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+(n - 2)Q VC 2 + (2 - n)2fo%C 2 - + 2(Q 2 - t^d 

' ' ((L t - Lq2 - 1)Q 4 + (Lq2 - L t )QH + t 2 ) + (L t - 4)t 2 



(t-Q 2 ) 2 \co 

+ (4 - 3L Q2 - Lq 2 + 2L t + L 2 )Q 4 + (L 2 Q2 - 3L Q 2 + 3L t - L*)QH 

tyl t>2 (2 ( 4 2 2 \ 2 



+ (t _ g2)3 " 1 " V)Q 4 + (L Q2 - L t )QH + t 2 ) + (2L t - 5)t 



+ (5 - 5L Q2 - Lq2 + 3L 4 + L 2 )Q 4 + (Lg 2 - 5L Q2 + 5L t - Lf)QH 

+ (148) 

where we have used the renormalization conditions for -Bq(O). Vertices of diagram V4 
and Vg are in this category. 

Here is a simple application of these methods in the case of a qqi * (at qx — 0). An 
explicit calculation gives to lowest order (in d = 4 — 2e dimensions) (with x = Q 2 /s) 

47r 2 o:e 2 . 

a° m ,( Sl , s 2 ) = -—^-8(1 - x) [(1 - /i) - e(/i - 1) - es 1± ■ s 2± ] . (149) 

oS 

Notice that, when we send Si ± and s 2± to zero, we reobtain the usual longitudinal 
helicity dependence of the matrix element. We observe that this cross section, away 
from d = 4, does not conserve helicity and has a spurious dependence from transverse 
spin, due to the e-dependent terms in ( |149| ). 



This violation originates from the t'Hooft-Veltman. Notice also that the cross sec- 
tion shows no dependence on the transverse spin at d — 4. Therefore, any transverse 
spin scatters, while the longitudinal components of the spin four-vector have to satisfy 
an helicity selection rule (the (1-h) factor). 

As in the longitudinal case, away from 4 dimension the t'Hooft Veltman prescrip- 
tion seems to indicate that there is a transverse spin dependence to lowest order, 
which is clearly an artifact of the regularization. 

Using the methods of this appendix, one gets quite straightforwardly for the virtual 
corrections to the (T qqi * cross section 

v An 2 aeQ 
a qqr (s 1 ,s 2 ) = — — -6(1 - z) [(1 - h) - e(h - 1) - esi ± ■ s 2 _ 



x ^ CF f^y_i_(-4-i- 8+ ^). (150 ) 



3s 

2n~ r \ Q 2 J r[l-e] V e 2 e 
19.4 gluon-quark-quark vertex (type 2) 

yn (27T/x) 4 -" r dni r(/+hh x ((i + 2p 2 y g ^ + (i - P2 ) y* - (21 + P2 y g ^) 

m 2 J / 2 (/+pi) 2 (/ + p 2 ) 2 
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We get 

V» = 2B^+i> 2 ^{2C l + C 2 + 2C,)-rhi>2{Ci-C 2 + C,)+YplC2 

+(n - 2) J* 2 (p£(C 2 + 2C 22 ) + 2rfC 12 ) + 2{n - 2)^C 00 . (152) 

1) In the case p\ — 0, p\ — Q 2 and (pi — p 2 ) 2 = we get 

V» = 2B Y + 2Q 2 ^(C 1 + C 2 + C ) + 2(n - 2) 7 ^ 00 (153) 

where 

5 = S (0), (154) 
C = 0,(0, 0,Q 2 ). (155) 

The vertex correction is 

V^ = y (1-2 + £q 2 ). (156) 

2) Pi = pi = and t = (pi — p 2 ) 2 = p\. 



V» = 2B 1 Y-tY(2C 1 + C 2 + 2C )-2p1 fa(Ci-C 2 + C ) 

+(n - 2)fo(p£(C 2 + 2C 22 ) + 2rfC 12 ) + 2{n - 2) 7 ^C 00 , (157) 

where 

B x = B (t) 

C = 0,(0,0,*) (158) 

and 

C 2 = -2C 22 . (159) 
Substituting these expressions and the renormalization conditions for B (0) we get 

„/2 2L 4 -3 7T 2 „\ ^ 2 „ / 1 ^ 7T 2 L, 2 \ , . 

\u z lu 6 J t \ur oj 12 2 J 

This vertex correction appears in the diagrams V3 and V§ of Fig. [7]. 

Notice that in all of the expressions given above, all the logarithmic functions 
containing a functional dependence in Q 2 or s of the form ln(— Q 2 /fi 2 ) or ln(— s//j, 2 ) 
should be analytically continued as ln(Q 2 //z 2 ) + in and ln(s//x 2 ) + m respectively. 
Notice that due to cancelations of the imaginary parts because of interference, we 
usually get ln 2 (— Q 2 / /j 2 ) — > ln 2 (Q 2 //i 2 ) — vr 2 and similar ones with Q 2 replaced by s. 
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20 Box contributions and tensor reductions 



20.1 massless box diagram 



The diagram is evaluated in the Euclidean region (s < 0, t > 0) and analytically 
continued to the physical region using the typical relations 



| s |e e «7rO(s) 



ln(— s) — > In |s| — 27r0(s) 



(161) 



where Q(x) is the step function. There are various expressions of the massless scalar 
box, with/without Spence functions which we list below. Their equivalence can be 
easily shown by using the relations 



Li 2 {x) + Li 2 (l -x' 1 ) = --ln 2 (x), x > 
ln(l — x) — ln(x — 1) + in. 

d 4 ~ 2e p 

(2vr) 4 -Vb - h) 2 {p -h- k 2 ) 2 {p - k x - k 2 - k 3 ) 2 



(47r) 2 - e st 
i 



((-sy + (-ty)-ln 2 (s/t)-n 2 



-s\ 



where 



cr 



(4ir) 2 st \AttJ 

,(_ 47r )(4-n)/2 

r[i + e ]r 2 (i- e ) 

r(l-2e) 



T[l + e] 



4 2,/ 5 2 
— mt/s — -7T 



"4 - n 


5 


'n 


n 






- 1, - - 1 


2 




.2 


2 



2(n-3)Ai 



A x = - (V e + r e ) - s -(i + ~) j(i + £)- r e (i + -) j(i + - 

st \e t t s s 

J{x) = -^^ + e -U 2 {x). 

X X 



(162) 



(163) 



(164) 



20.2 massive box diagram (scalar) 

The expression of the direct scalar box diagram with one external mass has been given 
in |TJ (here n = 4 - 2e) 



h(s,t,Q 2 ) = -i{4Ti) 



-n/2_ 



'{-Q 



2\-e 



-S) , U . 



st 

(-*) 



-Ffl,-e,l-e, 



nQ 2 
si 



, F(i,-6,i- e ,--; 

st s 



(165) 
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This expression can be cast into a simpler form by some manipulations which we are 
going to discuss here briefly. We use the expansion 



F(l, -e, 1 - e, z) = 1 + e ln(l - z) + e 2 S(-z) 



S(z) 



z , ln(l + x) 
ax 



-Li 2 (-z) 



(166) 



S(z) is the Spence function, related to the Euler dilogarithm. Expanding in e Eq. ( |165|) 

we get 

I 4 (s, t, Q 2 ) = -i(A-K)- n ' 2 c T {% \(-Q 2 y e - M~ e - Hrl + 2 lo gl , (167) 



where 



log! = - ln(-Q 2 ) ln(l + -^-) + ln(-s) ln(l + -) + ln(-t) ln(l + 

si r s 



S = 5 



S - )-S - 



u 



(168) 



We use the relations 



7T 



Li 2 (z) + Li 2 (l - z) = - ln(l - z) \nz + — , 

6 



Li 2 (xy) = Li 2 (x) + Li 2 (y) - Li 2 
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to relate dilogs of different arguments in order to get 
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After analytic continuation we get 



log! + log 2 
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and the final expression 
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(173) 



Notice that this relation needs analytic continuation in order to be applied to the 
calculation of the scalar box diagram. Its range of validity is in the unphysical region 
Q 2 < 0, s < and t < 0, therefore we analytic continue it to the region Q 2 > 0, s > 
and t < with the relations 



Q 



2\-e 



Is IV" 



(174) 



and leave t unchanged. Notice that there are no imaginary parts generated by this 
continuation, as one can easily check from ( |1 73| ) . 

In the case of the scalar w-channel diagram Vu we replace s — > u in (|173 ) and 
therefore we analytically continue only in Q 2 . There are imaginary parts generated 
by this procedure, which in this case don't cancel. However, since the interference 
terms are of the form Ab ox x B* Born + A* hox x Bsom in the differential cross section, 
than we get again a cancelation of the imaginary parts. The structure of the di-logs 
contribution generated by ( |1 73|) is different from ref. [[|. We have used 
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(175) 



to further simplify the result for the virtual corrections. 
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Figures 




Figure 4: Diagrams which contribute to the process g + g-^7* + g + g 
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Figure 5: Diagrams which contribute to the process q + q — > 7* + G + G 
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Figure 7: Diagrams which contribute to the process q + q — > 7* + G 
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